Transverse Fluctuations in the Driven Lattice Gas 
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We define a transverse correlation length suitable to discuss the finite-size-scaling behavior of 
an out-of-equilibrium lattice gas, whose correlation functions decay algebraically with the distance. 
By numerical simulations we verify that this definition has a good infinite-volume limit indepen- 
dent of the lattice geometry. We study the transverse fluctuations as they can select the correct 
field-theoretical description. By means of a careful finite-size scaling analysis, without tunable pa- 
rameters, we show that they are Gaussian, in agreement with the predictions of the model proposed 
by Janssen, Schmittmann, Leung, and Cardy. 



While the statistical mechanics of systems in thermal 
equilibrium is well established, no sound framework is 
available for nonequilibrium systems, although some in- 
teresting results have been recently obtained jTj. At 
present, most of the research focuses on very specific 
models. Among them, the driven lattice gas (DLG), in- 
troduced by Katz, Lebowitz, and Spohn §), has attracted 
much attention since it is one of the simplest nontrivial 
models with a nonequilibrium steady state || . The DLG 
is a generalization of the lattice gas. One considers a 
hypercubic lattice and for each site x introduces an oc- 
cupation variable n x , which can be either zero (empty 
site) or one (occupied site). Then, one introduces an 
external field E along a lattice direction and a gener- 
alization of the Kawasaki dynamics for the lattice gas 
with nearest-neighbor interactions. In practice, one ran- 
domly chooses a lattice link (xy), and, if n x ^ n y , pro- 
poses a particle jump which is accepted with Metropolis 
probabilities w(/3AH + /3E£), where I = (1,0,-1) for 
jumps (along, transverse, opposite) to the field direction, 
w(x) — min (1, e~ x ), and AH is the variation of the stan- 
dard lattice-gas nearest-neighbor interaction 

H = -A^n x n y . (1) 

As usual, the parameter (3 plays the role of an inverse 
temperature. A nontrivial dynamics is obtained by con- 
sidering periodic boundary conditions in the field direc- 
tion. Indeed, in this case a particle current sets in, giv- 
ing rise to a nonequilibrium stationary state that is non- 
Gibbsian. 



At half filling, the DLG undergoes a second-order 
phase transition. Indeed, at high temperatures the 
steady state is disordered while at low temperatures the 
system orders: the particles condense forming a strip par- 
allel to the field direction. The two temperature regions 
are separated by a phase transition occuring at the crit- 
ical value Pc(E) depending on the field E. f3 c (E) con- 
verges to the Ising critical value (3i for E — > 0, and, in- 
terestingly enough, increases with E, and saturates at a 
finite value /3 c (oo) when E diverges. Such a transition is 
different in nature from the order/disorder one occuring 
in the lattice gas. For instance, it is strongly anisotropic, 
i.e. fluctuations in density correlation functions behave 
in a qualitatively different way depending on whether one 
considers points that belong to lines that are parallel or 
orthogonal to the field E. 

Some time ago Janssen, Schmittmann, Leung, and 
Cardy (JSLC) proposed a Langevin equation for the 
coarse-grained density field, which incorporates the main 
features of the DLG, i.e. a conserved dynamics and 
the anisotropy induced by the external field and should 
therefore describe the critical behavior of the DLG. By 
means of a standard Renormalization-Group (RG) anal- 
ysis, critical exponents were exactly computed in generic 
dimension d for 2 < d < d c — 5. This model predicts a 
strongly anisotropic behavior with correlations that in- 
crease with different exponents V\\ and v±_ in the direc- 
tions parallel and orthogonal to the external field . The 
corresponding anisotropy exponent A = [vn/v±) — 1 can 



be exactly computed finding A 



G0/3. Another re- 
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markable property of the JSLC model is that transverse 
fluctuations become Gaussian in the critical limit. Ex- 
tensive numerical simulations [^) have confirmed many 
predictions of the JSLC model, although some notable 
discrepancies still remain. The conclusions of JSLC have 
been recently questioned by Garrido, de los Santos, and 
Muiioz 0. They argued that the DLG for E = oo is 
not described by the JSLC model but is rather in the 
same universality class as the Randomly Driven Lattice 
Gas (RDLG). This model has upper critical dimension 
d c — 3. The RG analysis done in Ref. pj] leads to se- 
ries expansions for the critical exponents in powers of 
e = d c — d. In particular, A = 1 + 0(e 2 ), consistent with 
A 1 in d = 2 (the case we will consider in our numeri- 
cal simulations), quite different from the JSLC prediction 
A = 2. A second notable difference is that in the RDLG 
transverse critical fluctuations are not Gaussian. Appar- 
ently, numerical simulations of the DLG || are also in 
agreement with the RDLG scenario and consistent with 
A « 1 in two dimensions. 

In view of these contradictory recent results, new nu- 
merical investigations are necessary, in order to decide 
which of these two models really describes the DLG uni- 
versality class. In this paper, we will focus on the trans- 
verse fluctuations for the very simple reason that in the 
JSLC model there is a very strong prediction: transverse 
correlation functions are Gaussian so that, beside criti- 
cal exponents, one can also compute exactly the scaling 
functions of several observables. Therefore, it is possible 
to make stronger tests of the JSLC model. 

In order to test the JSLC predictions, we will inves- 
tigate the Finite-Size Scaling (FSS) behavior of several 
observables. For an isotropic model on a lattice L d , the 
FSS limit corresponds to t = l—(3/(3 c — > 0, L — > oo, keep- 
ing tL 1 / 1 * constant. This has to be modified for strongly 
anisotropic models [[| like the DLG. It has been argued 



that, for a geometry L\\ x L 



both combinations tin 1 and tL 1 / v± 



one has to keep fixed 
and therefore also 



the so-called aspect ratio Sa 



-1/(1+A) 



/L±. Then, if 



an observable O diverges at criticality as t z ° , in a finite 
lattice one has 



0(0; Lii , L) « L z °/^f (t-^ /L, S A ), 



(2) 



where we have neglected subleading scaling corrections. 
To simplify the notation we write L instead of L\_ and, 
below, £l for the transverse correlation length. In many 
numerical studies Eq. (||) has been tested for several ob- 
servables, but this can be a very weak test since several 
parameters, /3 C , Vx, Zq, and A, must be tuned in order to 
fit the numerical data. A stronger FSS test can be per- 
formed if one uses a suitably defined correlation length. 
In this case, using a transverse (infinite volume) correla- 
tion length , Eq. (§) may be written in the form 



In this equation j3 c does not explicitly appear and zq 
and enter only through their ratio, cancelling when 
one looks directly at the correlation length. One can also 
eliminate this unknown by considering the ratio of the 
observables on two different lattices with transverse sizes 
L and aL. In this case we have 



0(/3;a 1+A L lb aL) 
0(J3;L h L) 



f W;L h L) 

to I ^ ,t> A ,Ci 



(4) 



where only A has to be fixed a priori. 

In order to test Eqs. (||) and (|J), one has to define 
a finite-size correlation length. In the DLG this is not 
obvious because correlation functions always decay alge- 
braically at large distances JM. A parallel correlation 
length was defined in Ref. [ |ll| , but it suffers from many 
ambiguities (see the discussion in Ref. pj). The def- 
inition of a transverse correlation length is even more 
difficult, because of the presence of negative correlations 
at large distances ||, || . 

In this paper we propose a new definition that gen- 
eralizes the second-moment correlation length used in 
equilibrium systems. The basic observation is that the 
infinite-volume transverse the two-point wall-wall corre- 
lation function decays exponentially, so that a transverse 
correlation length can be naturally defined in the ther- 
modynamic limit. The extension of this definition to fi- 
nite volumes requires some care because of the conserved 
dynamics, which makes the two-point function vanish at 
zero momentum. Here, we will use the results of Ref. J12) . 
Given the Fourier transform G(q) of the two-point corre- 
lation function {n x no), we focus on the transverse corre- 
lation G±(q) = G({q\\ = 0, q± = q}) and define 
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U±fe) 



(5) 



where q n = 2 sin (to / L) is the lattice momentum. If the 
infinite-volume transverse wall-wall correlation function 
decays exponentially, G±(q) has a regular expansion in 
powers of q 2 and 



G- A 1 {q)^ X ^ + b 2 q 2 + 0(q 4 



(6) 



O (/3; L I, , L) w L'° ^ fo (Coo (/?) /L , S A ) . 



(3) 



where the coefficient b of q 2 naturally defines a transvers 
correlation length ^ x . We expect Eq. (|^) to hold also in a 
finite box. Then, starting from Eq. (|^) , it is easy to show 
that converges to £oo as L — -» 00, justifying our defi- 
nition of finite-volume correlation length. In the subse- 
quent analysis we consider £13 as the finite- volume (trans- 
verse) correlation length } L (T) = £(T; S A +A L 1+A 1 L). 
As in previous studies, we also define a finite-volume 
transverse susceptibility as xl = G±(2w/L). 

As we already said, in the JSLC model transverse fluc- 
tuations are Gaussian. This allows us to compute the 
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scaling function appearing in Eq. (|j) for £l and XL- If 
G±(q) is Gaussian we have 



F ( (z,S 2 ,a)= [l-(l-a- 2 )(2^)V] 1/2 , (7) 

with z = £l/L and F x (z, S2, a) — Fg(z, S2, a) 2 . We can 
also compute f%(x, S2), see Eq. (||), obtaining 



/ ? (x^ 2 ) = (4^ + l/x 2 )- 1/2 . 



(8) 



In this paper we want to make a high-precision test of 
the theoretical predictions of the JSLC model, Eqs. (0) 
and (||). We work in two dimensions at infinite driving 
field. Since we wish to test the JSLC predictions we fix 
A = 2 and consider lattice sizes with S2 ~ 0.200. The 
largest lattice corresponds to Lu = 884, L — 48. For each 
lattice size, we compute xl and £l for several values of 
(3 lying between 0.28 and 0.312. 

As a preliminary test, we verify that our definition of 
£l has a good thermodynamic limit. For this purpose, 
we introduce the following quantity (the reason for this 
definition will be presented below) 



T L {f3)=^\f3) 



(9) 



In Fig. [I] we plot tx(/3) versus 1/L 2 at several inverse 
temperatures (3. For each (3, tl((3) converges to a finite 
constant, showing that our definition has a finite infinite- 
volume limit. Moreover, the same result is obtained by 
using sequences of lattices with S 2 or Si fixed: the result 
does not depend on the way in which L|| and L go to infin- 
ity. As expected, when the temperature approaches the 
critical value, it is necessary to use larger and larger lat- 
tices to see the convergence to the infinite- volume limit. 
At [3 fixed we expect the convergence to become even- 
tually exponential in L. However, for lattices with S2 
fixed we observed an intermediate region of values of L 
in which t l is apparently constant. Such a region widens 
as (3 approaches the critical point and is therefore in ex- 
cellent agreement with the relation 



C 2 (/3)«t l (/3)=£Z 2 (/3) 



Att 2 L~ 2 



(10) 



in the FSS limit L — > 00, (3 — > (3 C , neglecting sublead- 
ing corrections that will depend in general on the cho- 
sen value of 52 and are particularly small for our choice 
S 2 » 0.200. Eq. @ immediately gives Eq. (|). Thus, 
the results presented on Fig. [3] are perfectly consistent 
with A = 2 and the JSLC prediction for f^(x, S2). 

We want now to make a precise test of Eq. (ffl). In 
Fig. U we report the results of our simulations for the 
ratio !;2l/£,l as a function of for lattice sizes with 

fixed S 2 , and we compare them with Eq. (H). We stress 
that the theoretical curve is not a fit to the data: there 
is no free parameter to be chosen! Even if the agreement 
is not perfect, we notice that the points closer to the 
theoretical curve correspond to larger lattices. 
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FIG. 1: tl for different inverse temperatures /3. Filled (re- 
spectively empty) points refer to geometries with aspect ratio 
S 2 (respectively Si) fixed. Here & w 0.200, Si w 0.106. 
Errors are smaller than the size of the points. 



Using the universal function F^(z, 52,2), we can ex- 
trapolate our data to infinite volume using the gen- 
eral strategy of Ref. [|l3). Correspondingly, we obtain 
(3 C = 0.31256(9) and verify that, for small t, ^ ~ t' 1 ' 2 
as predicted by JSLC. 
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FIG. 2: FSS curve for the transverse correlation length at 
fixed & ~ 0.200. Different symbols correspond to different 
lattice sizes: L = 16(o), 18(D), 20(A), 22(0), 24 (*L The 
solid curve is the function F^(z, S2, 2) denned in Eq. (M). 
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FIG. 3: FSS curve for the transverse susceptibility at fixed 
S 2 » 0.200. We use the same symbols of Fig. |. The 
solid curve is the function F x (z, S2, 2) = F$(z, S2, 2) , where 
F${z, S 3 , 2) is defined in Eq. (p 



We can perform the same test for the susceptibility. 
In Fig. |^ we report our numerical results for xl together 
with the theoretical JSLC prediction. We observe a good 
agreement between theory and Monte Carlo results. 

We have also measured the transverse Binder cumulant 
<7l, see Ref. Q for the definition. At the critical point 
we observe that <?l(/3 c ) ~ L~ 0/2 for L — + oo. Thus, the 
Binder cumulant vanishes at the critical point, again in 
agreement with the idea that transverse fluctuations are 
Gaussian. 

Further details of our analysis are presented else- 
where [ p"4[ . 

In conclusion, we have shown by means of a new FSS 
analysis and new Monte Carlo data that the critical 
behavior of transverse fluctuations in the DLG model 
is Gaussian. Indeed, the numerical data are in per- 



fect agreement with Gaussian FSS functions. Note that 
our results are stronger than those presented in previ- 
ous analyses. First, we check not only critical exponents 
but a full scaling function; second, in our comparisons 
there are no free parameters that can be tuned. Our re- 
sults are in perfect agreement with the predictions of the 
JSLC model. 
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acknowledge the support of the INFM trough the Ad- 
vanced Parallel Computing Project at CINECA. 
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